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TECBNICAL NOTE 3158 


A SUBSTETUTE-STRINaER APEROACH FOR INCLUDING 
SHEAR-LAG EFFECTS IN BOX-BEAM VIBRATIONS 
By William W. Davenport and Edwin T. Kniszewakl 


SIMMARY 


The use of the substitute-stringer approach for including shear 
in the calculation of transverse modes and frequencies of hox beams is 
discussed. Various thin-walled hollow rectangular beams of uniform wall 
t hi c kne ss are idealized by means of the substitute -stringer approach 
the resulting frequencies of the idealized structures are compared with 
those of the original beams. The results indicate how the substitute- 
stringer idealization could be made in order to yield accurate represen- 
tation of the shear-lag effect in dynamic anclysis. 


INTRODUCTION 


In dete rminin g analytically the natiiral transverse modes and frequen- 
cies of box beams, the influence of shear-lag effects may be of consider- 
able in 5 )oii;ance, as indicated by investigations such as those presented 
in references 1 and 2. An appealing solution to the problem of 
shear-lag effects in a dynamic analysis of a built-iq? box beam sucdi as 
that shown in figvtre l(a) wotdd be to idealize the box beam into a simpler 
structure which involves fewer components but has essentially the same 
shear-lag properties. The simplest such idealized structure is the well- 
known substitvrfce-stringer structure. 

The substitute-stringer idealization is used by Kuhn and Peterson 
(ref. 3) in static problems for obtaining the maximum stresses of shell 
structures. There is, however, no indication that the idealized struc- 
tures which have been defined for static problems woxild be effective in 
determining natural modes and frequencies. It is true that Anderson and 
Houbolt xised the siibstltute-stringer idealization (ref. 2) to account for 
shear-lag effects on the natural bending frequencies of box beams. Their 
primary purpose, however, was to demonstrate the magnitvide of the shear- 
lag effects; no investigation of the accuracy of the approach was 
presented. 
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The pirrpose of the present investigation is to indicate how the 
substitute-stringer idealization can be made (or, more precisely, where 
to locate the substitute stringers) in order that the dynamic behavior 
of the prototype and of the Idealized structure will be essentially the 
same. This purpose is achieved by conparing the bending frequencies of 
several thin-walled rectangular tubes (which are analyzed exactly in 
ref. 1) with the frequencies obtained by an exact analysis of their ide- 
alized structures . 

In this paper, the idealization of an actual box beam into its 
substitute-stringer structure is discussed. The aforementioned compari- 
sons are then made and, conclusions are drawn with regard to the accuracy 
of the procedure. A list of symbols is contained in appendix A and a 
vibration analysis of the siibstitute-strlnger structure is included in 
appendix B. A pertinent extension of the solution of reference 1 is 
made in appendix C. 


THE SUBSTITUTE-STRIHGER IDEALIZATION 


A box beam which is typical of aircraft construction and its 
siibstitute-stringer stmcture are shown in figure 1 . The idealized struc- 
ture consists of four flanges and four stringers which carry only normal 
stress connected by sheets Tdiich carry only shear. The cross-sectional 
areas of the flanges and stringers of the idealized structure are deter- 
mined so that their moments of Inertia are the same as the moments of 
inertia of the spars and covers, respectively, of the original structtire; 
the moments of inertia in each case are taken about the horizontal axis 
of symmetry. The over-al 1 . dimensions and the web and cover-sheet thick- 
nesses are the same for both stnctures. The chordwise location of the 
substitute stringers, given by bs, is, however, as yet unspecified/ the 

value of bs determines the magnitude of the shear-lag effect in the 

idealized structure and is the quantity of paramount interest in this 
paper. Hereinafter, attention is directed to the effect on the vibration 
frequencies of varying bs and to the selection of the value of bg which 

yields accurate results. 


IX)CATION OF THE SUBSTITUTE STRINGERS 


In this section, conparisons are between the frequencies of 

variovis thinrwalled rectangular tubes such as that shown in figure 2 and 
the frequencies of their corresponding idealized structures. In each 
case the value of bg is pennitted to vary between 0 and b. The 
frequencies of the rectangular tiibes are obtained from a modification of 
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the exact series solution in reference the frequencies for their ide- 
alized structures are obtained from the solution presented in appendix B. 
The exact series solution of reference 1 is modified to include only the 
secondary effects considered in the substitute-stringer solution of appen- 
dix B, that is, shear lag and transverse shear deformation. 

The web and cover-sheet thicknesses and the cross-sectional areas 
of the flanges and_ stringers of the idealized structure are obtained as 
outlined in the preceding section and are (see figs. l(b) and 2) : 

tw = t(j = t 
Aj. = i at 

Al = bt 

In 02Tder to preserve the inertial properties, the mass per unit length 
of the substitute-stringer structTire is taken eq.ual to that of the rec- 
tangular tube. 

The effect of varying bg on the frequencies of the first three 

symmetrical free-free modes of the idealized structTire is shown graph- 
ically in figures 3 to 6. In each case, the frequency cd is expressed 
in the form of its relative error - 1 vhen compared with the exact 

frequency of the rectangular tiibe, and b3 is expressed in the form 

of the ratio bg^b^, wherein bQ is the distance from the web to the 
centroid of area of the half-cover. This ratio bg^bp is used to accord 

with past practice in static shear-lag investigations and also in the hope 
that the results will be applicable to more general types of box beams . 

It should be noted that, for the rectangular tiibes, bg equals b/2. 

The first case considered is that of a rectangular tube with a 
cross-sectional aspect ratio b/a of 3»8 and a plan-form aspect ratio 
L/b of 6.0. The curves of flgtire 3 cross the line of zero error at dif- 
ferent value's of tg^b^; thus no single value of bg^bj, gives exact fre- 
quencies for »n the modes. It is possible, however, to choose an "opti- 
mum" value which gives nearly exaet results for all the modes considered. 
This optimum value, ’irtiich has arbitrarily been selected so that the maxi- 
mum of the errors in the freqviencies of the first three modes is a mini- 
mum, is bg^bg = 0.56. The ^nfl.v^■m1Tm error in the frequencies for this 

value is less than 1 percent. 
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The effects of different cross-sectional and, plan-form Eispect ratios 
are indicated, in figures li- and 5* In figure 4, results are shown for 
cross-sectional aspect ratios of 1 and » (a limiting case) with a plan- 
form aspect ratio of 6.0. In figure 5 , results are shown for plan-form 
aspect ratios of 2.0 and 12.0 with a cross-sectional aspect ratio of 3«6. 
The CTurves in figures 4 and 5 are similar in character to those presented 
in figure 3; they differ only in steepness and in the values of tg/bg 

where the zero crossings occur. The optimum VEO-oae of and the 

maximum percentage error of this value for each case are included in 
table 1. For all these cases, the maximum percentage errors are very 
small. It should he noted that, except for the case in which b/a = 1.0 
and L/b = 6.0 and that in which h/a = 3*8 and L/b = 2.0, the optimum 
values of bg/b^ fall within a small range. In the first of these two 

cases, the shear-lag effect is very small in magnitude (see fig. 4); the 
second case is an extreme configuration, having a plan-form aspect ratio 
of 2.0. The results indicate, therefore, that for reasonable configura- 
tions with appreciable shear-lag effect, the optimum value of bs^bg is 

relatively independent of the cross-sectional and plan-form aspect ratios. 


The rectangular tubes treated thus far are admittedly not very real- 
istic, and the extensibility of the results obtained to more usual box 
beams, such as that shown in figure l(a), is questionable. For this 
reason, a generalization of the rectangular tube which more nearly rep- 
resents actiial structures has been considered. This generalized rectan- 
gular tube is assumed to have at each point of its cross section a thick- 
ness t that carries shear and also a different thickness t' that 
carries no rm al stress; this assimrption approximates the sitviation in an 
act\ial structvcre in which flanges and stringers carry normal stress but 
do not carry shear. The exact series solution of reference 1 can be 
extended to this "dual-thickness” staructure by the modification of the 
parameters shown in appendix C. 


A particular exaniple of this type of structure is considered in 
figure 6, where L/b = 6.0, b/a = 3*6, and t' ft = 2.0. Once again the 
mflvinnTm error for the optimum value of bg^bg is small (see table l) . 

The optimum value of Is less than the value obtained for the 

single-thickness comterpart; this reduction indicates a possible depend- 
ence of the optimum value of on the magnitude of the additional 

stress-carrying area of an actual structure. 
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EVALUA330N OF RESULTS 


The ratios lised for the configurations presented herein bracket 
those 'which 'would occur in most actual box beams. Cross-sectional and 
plan-form aspect ratios have been 'varied from 1 to oa and from 2 to 12, 
respecti'vely. The ratio t’/'t^ 'which is a measure of the ratio of normal- 
stress-carrying area to shear-carrying area in the co-ver, has been 
'Varied from 1 "to 2^ representati've 'values of t‘/i^ for actual box beams 
are 'within this range. The results in table 1 Indicate that althou^ 
these -wide ranges of ratios -were used, the optimum -value of "bg^bg falls 

-within a relati-vely narrow range. This result suggests -fche possibility 
that, in the absence of a more detailed investigation, a uni-versal -value 
of b^bg could be used. In order to demonstrate the -validity of -this 

conclusion, -the -values of the maximum of -fche errors in the first three 
-vibration frequencies for a -value of bg^bg = l/2 are eOso presented 

in table 1. The errors exceed 2 percent only in the extreme cases. 


Langley AeronaTifcical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., No-vanber 2, 1953 • 
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APPENDIX A 


SmBOLS 


cross-sectional area of flange of stibstitute-stringer 
structure 

cross-sectional area of substitute stringer 
effective shear-carrying area 
parameter defined by equation (C 7 ) 


half -depth of beam, 
half -width of beam 

distance between web and adjacent substitute stringer 
distance between web and centroid of area of half-cover 


constant 

modiilus of elasticity 

shear modulus of elasticity (taken equal to E/2.65 in 
present paper) 

bending moment of inertia 

, / Gtn An} 

shear-lag parameter, \l — ^ 

Y Ebg AjAl 

parameter defined by equation (C 5 ) 


frequency coefficient, oj 


coefficient of shear rigidity, — \/ 

lVgas 
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T 


L 
Nn 
N ' 
Pn 

P 

s 

t 

■tc 


t’ 


T 

U 

u(x,s) 

w 

X 


half-length of free-free heam 
parameter defined hy equation (B22) 

parameter defined hy equation (C^t) 

parameter defined hy equation (B23) 

perimeter of cross section of heam 

distance along periphery of cross section of heam 

wall thickness of heam with uniform wall thickness; effec- 
tive shear thickness for dual-thickness heam 

cover-sheet thickness 
weh thickness 

effective thickness for normal stress for dual-thickness 
heam 

maximum kinetic energy 
maximum strain energy 

longitudinal displacement of a point on heam 

longitudinal displacement of point of flange in siibstitute- 
stringer structure 

longitudinal, displacement of point of stibstitute stringer 

vertical displacement of cross section of heam 
longitudinal . coordinate 
Fourier series coefficients 


l,n 




P2 


integers 

1 

> parameters defined hy equations (B 26 ), (B27), (B29)^ 
and (B30) 
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e 

V- 

(0 




°L 

TC 


•^w 


incUiiation of normal to the wall with the vertical 

mass of heam per unit length 

natural frequency of heam obtained hy use of substitute- 
stringer approach 

natural frequency of heam obtained by use of exact solu- 
tion of refer^ce 1 

longitudinal stress in flange of substitute-stringer 
stmcture 

longitudinal stress in substitute stringer 

shear stress in cover sheet of substitute-stringer 
struetTrce 

shear stress in web of substitute-stringer structure 



APPENDIX B 


VIBEIATION SOLUEON OF A FREE-FREE 
SUBSTITUTE-STRINGER STRUCTURE 


The natural inodes and frequencies of the siibst.itute structiire may 
he obtained by the method employed in reference 1, that is, the Raylei^- 
Bitz energy procedure in conjunction vith appropriate Fourier series 
expressions . 

Let X be the longitudinal coordinate with its origin at the mid- 
point of the beam; then, by Hooke's law and the assumptions concerning 
the stress-carrying properties of the components of the structure given 
in the body of the paper, along with the assumption that cross sections 
maintain their shapes, the longitudinal stresses are 


ffp = 


dup 

dx 


(Bl) 


dUr 

Ot — E 
^ dx 


and the shearing stresses are 


rc = G 


= G ^ 

” Vdx a J 


(B2) 


(B5) 

(B4) 


where vtp is the longitudinal displacement of a point on a flange, 

is the longitudinal displacement of a point on a stringer, and w is 
the vertical displacement of a cross section of the beam. 


From these expressions for the stresses, the maximum strain energy 
of the structure is 
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and the maxinrum kinetic energy Is 

pL 

T = — (b6) 

2J-L 

where and w are now considered as the anrplitudes of displace- 

ment for the particular mode considered, a> is the natural frequency 
of the mode, and ji is the mass per unit length of the structure -vSiich 
the substitute structure represents. 

A natural mode of vibration must satisfy the variational equation 

5(U - T) = 0 (B7) 


where the variation is taken Independently with respect to Uj,, and 

w. Application of this principle to expressions (B5) and (b 6) would result 
in the differential equations anfl the natural boundaiTr conditions of the 
vibrational problem under consideration. However, Fourier series expres- 
sions for t^, and w are used in conjunction with the variational 

procedure, rather than a direct attack on the differential equations and 
the boundary conditions. 

Appropriate assunptions for the displacements for the symmetrical 
modes of a free-free beam are 


00 

w = C + % cos ^ 

n=l,5,5 ^ 

(B8) 

^F = £ ^n sin 

n=l,5,5 ^ 

(B9) 

^ = £ °n sin ^ 

n=l,5,5 2 L 

(BIO) 

Sxibstituting the expressions (b 8), (B9), and (BIO) into equations (B5) 
and (b 6) and then \ising expression (B7)^ where the variation is with 
respect to the a's, b's, c’s, and C independently, results in the 
following equations (where i = 1, 5^ 5# • • •) : 
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iff - J f - 1=1 - 21^ V ^ (- 1 ) 


1-1 

2 


C = 


ftfSlsta, + 

L Ap 2 ^ 


* ==s"(KL)" ^ * 



0 (Bll) 
= 0 (B12) 



|(XL)]c, 


0 



E ^=n(-l) 


n-1 

2 


= 0 


(BI 3 ) 

(Bl4) 


These equations are -written in terms of kg, the coefficient of shear 
rigidity, K, the shear-lag parameter, and kg, the frequency coefficient 
these terms are defined as follows: 


= hi— 

^ l]/ GAs 

(B 15 ) 

g- _ 1 Mtc At 
yEbs ApAi, 

(B 16 ) 

8 

II 

(B 17 ) 

II 

1 

(B18) 

I = 

(BI 9 ) 

At = Ap + Al 

(B20) 


By sol-vlng equations (BU), (B12), and (BI 3 ) for a^ and substituting 

the result in-to equation (B14), the following frequency equation, which 
must be satisfied by the frequency coefficient kg, is obtained: 
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i+^bV 2: (if)‘ 

n=l,3,5 ' ' 


- st)‘ 


= 0 


(B21) 


where 




(B22) 


W) * 


Pn-(f) -J^bV (“5) 

The rate of convergence of the series of equation (B21) is increased 

CO ^ 

by STibtractiiig the e3q)ression 2 ^^^^ ^ ^ adding the 

n=l,3,5 (^)^n 

p/tan hrjJEc! \ 

equivalent closed-form expression I s_2, _ ij . The restating 

\ ^b’^ / 


equation is 




tan + 2%^’ Z1 ^ 


= 0 


(B2lt) 


The series in equation (B21) converges as l/n^ while the series in 
equation (B2l^) converges as l/n^. 

For the case where h/a — >«, equation (B21) reduces to the closed 

form 




(B25) 
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where 


and 


ai 


1 Ebgh 

2 GL^ 


V 

2 



+ 



(B 26 ) 



(B27) 


For the case irtiere hg = 0, equation (B21) reduces to the closed form 


where 



tanh kgot^. + ^B^2^ " ® 


(B28) 



This frequency equation is identically the equation obtained when 
the frequency equation for symmetrically vibrating free-free beams of 
reference 4 is modified by neglecting rotary inertia, that is, when the 
only secondary effect considered is transverse shear deformation. 

Mode shapes of the structvcre may be obtained by solving equations (BU), 
(BI 2 ), and (BI 5 ) for a^, b^, and and then substituting the results 

into escpressions (b 8 ) , (B9) j and (BIO ) } the value of C may be arbitrar- 
ily chosen. 

Esctension to the antisymmetrical modes and to cantilevered beams 
may be accomplished by methods similar to those shown in reference 1 . 
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APPENDIX C 

EXACT FREQUENCY SOLUTION POE A DUAL-THICKNESS EECTANGULAE TUBE 


The exact solution of reference 1 can he extended to tape Into 
account the effects of having different normal-stress- and shear-carrying 
thicknesses t' and t hy modifying the expression for the maximum 
strain energy as follows : 



The expression for maximum kinetic energy is unchanged if longitudinal 
inertia is neglected and is 


T 



(C2) 


By means of the pirocedure described in reference 1 , the fol3x>wing 
frequency equation can he obtained for a symmetrically vibrating, free- 
free "dual-thickness" rectangular tiibe: 


^B 


1 + ki 


:.a 12 f-)® — 


= 0 


(C3) 


where 


N ' 
“n 


n^ 

8ks^ 


'jt - 


K* p 

P a 
Ifks n 


sum « - x) 

2 K'\P / 


cosh — ^ 
2 K' 


+ tanhSL^ 
2 K' 


ik.2 


(CM 
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and 

K' = il/II 

(C5) 


ho = 

^ L VgAs 

(C6) 


Ag* = l|-at' 

(C7) 


II 

(C8) 


It should he noted that, with the exception of the slight c2iange in 
parameters, equations (C3) and (C^l-) are the same as equations (^l) and 
(A18) in reference 1 if, of course, the effects of longitudinal inertia 
are neglected. 
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TABLE 1 

A COMPARISON OP ERRORS FOR OPTIMUM VALUES OF “bg/bc 
AND FOR A VALUE OF bg/Bg of 0.5 


Plan-form aspect ratio L/b = 6.0 

Cr os s -se ctlonal 
aspect ratio, h/a 

Optimum value 
of bg/bg 

Maximum percentage error for - 

Optimum value 
of bgyTbg 

^s/^c = 0*5 

1.0 

0.65 

0.5 

1 

* 3.6 

.56 

1 

2 

00 

.5^ 

2 

2.5 

Cross-sectional aspect ratio b/a =3-6 

Plan-form 
aspect ratio, l/b 

Optimum value 
of bg/bg 

Maximum percentage error for - 

Optimum valxie 
of bg/bg 

^s/^c “0*5 

2.0 

0.14 

2.5 

5.5 

*6.0 

.56 

1 

2 

12.0 

.61 

.5 

1.5 

Cross-sectional aspect ratio b/a = 3-6; 
plan-form aspect ratio L/b = 6.0 

Thlclmess ratio, 
t*/t 

Optimum value 
of bg/bc 

Maximum percentage error for - 

Optimum value 
of bg/bg 

^sA*c - 0*5 

*1.0 

0.56 

1 

2 

2.0 

.48 

1.5 

2 


This case is repeated for ease of comparison. 
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(a) Typical box beam. 



Rgure L- Typical box beam and its substitute-stringer idealization. 
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Rgure 3 .- Effect of stringer location on the accuracy of the substitute 
stringer approach for a box beam of uniform wall thickness. 



Rgure 4.- Effect of stringer location on the accuracy of the substitute-stringer approach for box beams with extreme 


cross-sectional aspect ratios. 
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be 


Figure 6. “ Effect of stringer location on the accuracy of the substitute- 
stringer approach for a dual-thickness box beam ( = 2). 
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